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Abstract 

Theory of the first order Raman scattering (FOR) line shape allowing for nonlinear 
and correlation effects contribution to inhomogeneous broadening as well as dynamic 
mechanisms of homogeneous broadening is developed. It is shown that in general case 
FOR line contains two maxima, shape of low frequency one being defined by homoge- 
neous broadening. Our theory explains the peculiarity of observed FOR scattering of 
TO2 hard phonon in KTL and KTN at different temperatures and concentrations of 
Li and Nb ions. 



1 Introduction 

Investigation of Raman spectra of the first order (FOR) in disordered ferroelectrics was 
shown to be sensitive method for studying of the critical slowing down of optic phonons 
nearby ferroelectric transition and dynamics of local fluctuation (see e.g. and ref. 
therein). Several attempts were made to explain FOR spectral anomalies such as strong 
increase of peak intensity in the vicinity of T c and peculiar shape of the line || . How- 
ever, up to now there is no general approach to the description of Raman spectra in the 
disordered ferroelectrics nearby T c . These results are in the discrepancy in the values of the 
materials physical characteristics obtained from experimental data (see e.g. flj] -@|). The 
most promising way for development of general approach seems to be the model allowing 
for inhomogeneous broadening of FOR lines induced by static disorder and homogeneous 
broadening due to dynamic effects. The existence of both these contributions to Raman 
spectra was shown early in KTa\- x Nb x O^ (KTN) and K\^ x Li x TaO^ (KLT) fl], |2|, || . 

Theory of inhomogeneously broadened lines in radio-, optic- and other spectroscopy meth- 
ods was developed in many details early for conventional dielectric and magnetic systems 
0, 0. In these materials inhomogeneous broadening used to be temperature independent, 
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so that only homogeneous broadening may be the source of temperature dependence of the 
spectra line shape and width. In the disordered ferroelectrics due to nonlinear and spatial 
correlation effects which is known to be especially strong in the vicinity of T c inhomogeneous 
broadening was recently shown to be temperature dependent ]F| . 

In the present work a theory of FOR allowing for both linear and nonlinear contributions 
to inhomogeneous broadening as well as to homogeneous mechanisms of line broadening is 
developed. The theory explains the main peculiarities of observed FOR spectra in KTL with 
1% and 4% of Li || and in KTN with 15.7% of Nb ions when the systems approach the 
phase transition from above. Since in vicinity of the phase transition the nonlinear effects are 
known to be large enough our theory is more general than that in [0, [3] , where nonlinear and 
spatial correlation effects were not taken into account. The proposed theoretical approach 
allows to describe of the Raman spectra (in distingvish 0, |3j ) without calculation of the 
correlation function of the quasi-static polarization fluctuations but only suggesting about 
its spectral density form. Because of generality of proposed theory it can be applied to many 
ferroelectrics in vicinity of transition temperature. 



2 Theory of first order Raman spectra shape 

2.1 Incipient ferroelectric KTaO^, doped by Li or Nb ions can have ferroelectric phase tran- 
sition, mixed ferroglass phase and dipole glass state in dependence on temperature and the 
impurities concentrations [§, |9|. These materials can be regarded as the model systems of 
the disordered ferroelectrics. Their peculiar properties are connected with off-center posi- 
tions of Li and Nb ions substituted for K + and Ta 5+ ions respectively. Some anomalies of 
Raman spectra can be explained also by these ions off-centrality. In particular quasistatic 
polarization fluctuations induced by off-center ions can lead to the reduction of KTaO^ cu- 
bic symmetry and to the appearance of first order Raman scattering (FOR) even above the 
transition temperature T c [[]], [3] . The most detail description of the appeareance of the single 
phonon Raman lines above T c can be found in QT0j| where it was shown how the slow finite- 



range precursor order, dominant near the phase transition, is responsible for the persistence 
of quasi-first-order hard-mode features above T c . 

As known, Raman scattering in perovskite ferroelectrics is caused by the changes in the 
oxygen electrtonic polarizability 5a(r, t) induced by optic vibration modes of the lattice. 
This changes can be written in the form 

6a(r,t) =P(r,t) • A-P(r,t) (1) 

where P(r,t) represents the space- and time-dependent polarization fluctuation and A is a 
fourth-rank tensor. Equation (1) usually describes second-order Raman scattering. However, 
if we write the fluctuation polarization as a sum over the components due to the polar hard 
modes P h and slow relaxing component P M we can explain the appearance of single-phonon 
Raman scattering above T c considering the cross terms pvp h i n the Eq.(l) The scattering 
intensity is given by the spatial and temperature Fourier components of the polarizability 
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correlation function 



I{uj) ~ < 5a(r,t)Sa(0,0) > q=0tUJ 

~ E / ^' < ^ *) pM (0> 0) > gV < P ft (r, t)P h (0, 0) >_ gW , (2) 

Because the first correlation function has a sharp maximum near uj' — and second one near 
/(cj) is sharply peaked at uj = Q q , i.e., near the hard single-phonon frequency. 
It is seen that (Q) can be represented as the convolution of the functions represented 
inhomogeneous and homogeneous mechanisms of broadening connected respectively with 
the first and the second correlators and so it can be written in the form 

I(uj) = I J(uj,uj')f(uj')duj' . (3) 



The detailed form of homogeneuos J(uj,uj') and inhomogeneous contribution f(uj') can be 
obtained by the following way. Thermal lattice vibrations usually lead to homogeneous 
broadening with simple Lorentzian shape, i.e. 



< P h (r,t)P h (0,0) >_i^~ — ^ ; - — t ( I.) 



n q > r 2 + (u-uj>- n a ,) 



which should be valid for the difference of frequencies {uj — uj') close to hard phonon frequency 
Q q t. On the other hand the correlation function of quasistatic polarization reads 

< P»(r, t)P^0, 0) >m=< P^r)P"(0) > q , it6{uj') . (5) 

Substituting (long-wave approximation) J2 q > — y j^js I d 3 q' an d performing some transforma- 
tion one obtaine that 

/(u/) = J dq'q' 2 < P»(r)P»(0) > q , 5{uj' - fy)^- , (6) 

r v 

1 1 + [uj — uj' 1 ) In 

Equations (^|) and (0) determine respectively the contribution of inhomogeneous and homo- 
geneous broadening into Raman line shape (see Eq. ([3])). However it is seen from Eq. (|B|), 
that for f(uj) calculation one has to calculate firstly the correlation function of polarization. 
Its calculation utilized some special models and supposition |], || . To our mind more gen- 
eral approach for f(uj) calculation can be statistical theory of specrtal line shape ||. This 
theory was sucsessfully applied for description of optic, radiospectrroscopy, 7-resonance etc. 
line shapes || |TT|. 

2.2 The shape of f(uj) was calculated in the statistical theory framework for the cases 
when the frequency shift of spectral line is linear || ^| and nonlinear []7|, |Tl| function of 
random fields. The latter case seems to be especially important for disordered ferroelectrics 
at T = T c ± AT, AT ~ (10 4- 20)K, where nonlinear and correlation effects is known to be 
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large enough. As a matter of fact the distribution function allowing for these effects can be 
expressed via linear one in the framework of the general theory of probability which makes 
it possible to write down the distribution of one random quantity via that of another if the 
relation between these quantities is known | 12| . For example in the simplest case when the 
random quantity x is a single-valued monotonous function of another random quantity h 



then the relation between their distribution functions g(h) and f(x) can be written as [12] 



g(h) = f(x(h)) 



dx{h) 



dh 



In general case when several different x values correspond to the same h(x) value the space 
of x should be devided into the regions where the function h(x) is monotonous. For the 
entire x-domain distribution function g(h) can be represented as a sum of terms like Eq.(8) 
|7], [12] . Since the distribution function in the case of linear contribution of the random 
field fi{oj) can be analytically calculated in statistical theory approach [g, |J we have to 
express via it the distribution function allowing for nonlinear random field contribution. Let 
us suppose that the shift Auj = uj — ujq = uj of the spectral line maximum position ujO due 
to random field contribution can be written as a power law (up to some m th power) of the 
random field uj' 

oj = uj' - a 2 uj' 2 - ... - a m uj' m . (9) 

Equation (9) makes it possible to express the distribution function of uj via that of uj' in the 
form @, 0, U 

d<p(uj, uj') 



?c=i 



(uj, uj') = UJ 



[UJ- 



UJ 



duj' 



/2 

OL 2 UJ — 



- OL m UJ 



/in 



where uj^ are the real roots of the algebraic equation 



[OJ.UJk 



0. 



(10) 
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It is seen, that m th order distribution function / m (co>) is expressed via that calculated in 
linear approximation when all nonlinear coefficients equal zero, i.e. a 2 = «3 = • • • = ot m = 0. 
Shape of fi(uj) calculated in the statistical theory approach was shown to be Gaussian, 
Lorentzian or Holtzmarkian in dependence on types of random field sources with parameters 
determined by the sources concentrations and characteristics @. [| . 

In the simplest case, when the main contribution is connected with the first nonlinear 
term in Eqs. (D, (02) (a 2 ^ 0, a 3 = . . . = a m = 0) Eqs.(10), (11) and (12) lead to the 
following form of the normalized second order distribution function 



J2(UJ) 



a/1 + 4a 2 uJ 



h 



a/1 + Aa 2 uJ — 1 
2^ 



fi 



a/1 + 4a 2 uj + 1 
2a 9 



(13) 



where G is the teta-function, so that f 2 7^ only in region uj c < uj < 00 (a 2 > 0) or 

—00 < uj < uj c (a 2 < 0), 

uj c = -(l/4a 2 ) (14) 
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is a critical frequency at which a divergency of f 2 {oj) appears. Therefore f 2 (u;) is strongly 
asymmetrical, and its form at uj ^> uj c (the wing of f 2 {uj)) is the following 



f 2 (u) -> 1/ VI + 4 «2^ • 



(15) 



Let us consider the homogeneous contribution in the Lorentzian form (see Eq.(7)). In this 
case the integration in Eq.fQ) is equivalent to substitution uj ± 1 for uj in the Eq. (|I3|) (1/r = Y 
is half width on the half height). If there are several mechanisms of homogeneous broadening 
with Lorentzian forms 1/r = YU ^-/ T ii where i numerates the mechanisms. 

The aforementioned procedure (or integration in Eq.(|3] )) with respect to Eq.(|l"3"|) in 
supposition that f\(uj) has Gaussian form, i.e. 



AM 



a; 2 
'2A 2 



27T 



(16) 



leads to the following shape of spectral line in the considered case 



h + - arctanrfcj 

1 7T v 



Aj27T(p(uj) 



exp 



S^uj) -2(1 + 2a 2 uj) 



+ exp 



where 



8a 2 2 A 2 

-S 1 (u)-2(l + 2a 2 uj) 
8a 2 2 A 2 

(1 +Aa 2 u) 2 + (4a 2 /r) 2 : 



cos 



4a 2 /r - #2(0;) 



cos ■ 



8a|A 2 

Aa 2 /r + S*2(u;) 
8a 2 A 2 



(17) 



S 1>2 {uj) = V2y/<p(u) ± (1 + 4a 2 cu) . 

The results of numerical calculations for several values of dimensionless parameters a 2 A and 
1/(tA) are depicted in figs. 1,2. It is seen that homogeneous contribution transforms f 2 {uj) 
divergence at to = u c into sharp maximum, so that the spectral line has two maxima (see 
fig. 1) instead of one in the linear case origin of high frequency one being connected with 
Gaussian form. The distance between two maxima approximately equals u c (see fig. 1). At 
large enough nonlinear contribution only sharp maximum at uj = uj c conserves, its width 
increases with 1/r increasing (see fig. 2). The left hand side half-width of low-frequency 
peak completely defined by homogeneous broadening mechanisms, meanwhile right hand 
side part of the line at uj > uo c defines mainly by inhomogeneous mechanism contribution. 



3 Raman spectra in KTL and KTN 

3.1 The developed theory was applied to recently observed T0 2 FOR in KTL with 1% and 
4% of Li and in KTN with 15.7% of Nb || 0. 

Let us begin with the consideration of Raman spectra in KTL. Measurements were 
carried out at T = 10 K (xu = 0.01) and T = 55 K [xu = 0.04) 0. In both samples 
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the observed line was strongly asymmetric with maximum at uo ~ 198 cm 1 . These lines 
were fitted good enough by Eqs. (12) - (14) with the following dimensionless parameters 
a 2 A = 0.38, (rA)- 1 = 0.06 (fig. 3) and a 2 A = 0.55, (rA)- 1 = 0.15 (fig. 4) for the considered 
Li concentrations respectively. 

The position of the lines maxima u> m — ujq = u c is defined by the parameter of nonlinearity 
in accordance with Eq. ([14]). We obtained a 2 — 0.1 cm by the fitting of the high frequency 
line tails with the Eq.(15). This value gives uj m = u c + u) Q ~ 197.6 cm -1 with the reference 
point u = 199 cm -1 . Note that Eq.(15) describes also observed frequency dependence of 
the lines tails for the same a 2 value. This made it possible to obtaine yhe magnitude of the 
Gaussian width A with the help of aforemationed a 2 A values: A = 3.8 cm -1 [x Li = 0.01) 
and A = 5.5 cm -1 (x Li = 0.04). 

In accordance with our theory the line form at uj < uj c defines completely by homogeneous 
mechanism contribution represented by Eq. (4) . It is seen that low-frequency half- width 1/r 
is connected with hard phonon life time T -1 . Keeping in mind the obtained values of 
(rA)" 1 = T/A and A one find T (T = 10 K) ~ 0.25 cm -1 , T (T = 55 K) ~ 0.8 cm" 1 . These 
date are in resonable agreement with ordinary values of hard phonon life times and with 
observed low-frequency half- width (see figs 3, 4 ). 

3.2 Now let us proceed to consideration of Raman spectra in KTN with 15.7% Nb which 
has the transition from cubic to tetragonal phase at T c = 138.6 K. The measurements were 
carried out at several temperatures in vicinity of T c : T = 160, 150, 146 and 142 K ||. 
At all the temperatures the lines with two maxima were observed, the first being nearby 
u)\ ~ 200 cm -1 and the second at uj 2 ~ 220 cm . The intensity of low frequency maxi- 
mum essentially increased with T lowering, meanwhile the intensity of high frequency one 
became very small at T = 142 K. Since nonlinear parameter a 2 increases with temperature 



approaching to T c JTT|, nonlinear effects has to be the largest at T = 142 K. Qualitatively 



transformation of spectra from two-peak line to one-sharp peak line at nonlinear parameter 
increasing is in agreement with theoretical overcasting (see figs.l, 2) It is seen also that 
increasing of low frequency maximum intensity with T lowering can be the result of 1/r 
decreasing. 

To be sure that nonlinear effects are really responsible for observed Raman spectra trans- 
formation we checked if the Eqs. (13-15) fitted experimental spectra. It was shown that 
Eq.([Lj|) fitted the wing (u > 230 cm" 1 ) of Raman line at T = 142 K for a 2 ~ 0.015 cm, 
which lead to u c ~ 20 cm -1 (see Eq. flT4]) ). This value fits pretty good the distance between 
two peaks of observed Raman spectra (see fig. 5), which speaks in favour of nonlinear ef- 
fects contribution. In fig. 5a we depicted line shape calculated with the help of Eq. ( |i~7|) for 
a 2 = 0.012 cm and A ~ 17 cm -1 . The later quantity was taken from measured high fre- 
quency maximum width. Note, that obtained A value made it possible to fit dimensionless 
u)j A scale with u scale in fig. 5 and later in fig. 6. 

The values of homogeneous broadening parameter 1 /r at different T were calculated in 
supposition that it defines mainly by reorientational frequency of elastic dipole connected 
with Nb. This frequency temperature dependence was measured early |13| and it was de- 
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scribed by Arrenius law 

- = -e W (-U/T) 
r t 



with U=200 K and l/r = 7 x 10 9 Hz. At T = 150 K Eq.flT3|) leads to the value (rA)" 1 
0.004, which was used in calculation of the line, depicted in fig. 5a. It is seen from fig. 5 
that the calculation and observed spectra look like one another. More detailed comparison 
of the theory and experiment was carried out for Raman line observed at T = 142 K (see 
fig. 6). Theoretical curve was drawn for aforementioned parameter = 0.015 cm, obtained 
from line wing behaviour and 1/r was calculated with the help of Eq . fll8|) . It is seen that 
theory fits pretty good observed Raman spectra. This gave evidence that measured Raman 
line asymmetry and rapid change in line shape when approaching the transition from above 
is really connected with nonlinear effects in KTN with 15.7% Nb. Theory overcasts strong 
increasing of low frequency peak at the region 142 K > T > T c because of «2 increasing and 
1/r decreasing in supposition that there is no another even small temperature independent 
contribution to 1/r. 

Comparing the obtained data for KTL and KTN one can see that because the reori- 
entation rate of Nb dipoles is much larger than that of Li, the homogeneous broadening of 
Raman line in KTL is defined by hard phonon dynamic whereas in KTN - by Nb elastic 
moment reorintation. Note that in 0, |3|] the reorientation of Nb electric dipole moment were 
supposed to be the origin of the homogeneous broadening. However the parameters of the 
electric dipoles orientations obtained from fitting of the theory with the experiment were 



strongly different from those obtained early in independent measurements [14 



4 Discussion 

The proposed theoretical description of Raman spectra shapes was perfomed without cal- 
culation of correlation function of quasy-static fluctuation of polarization. The comparison 
with the caculations of Raman spectra based on calculation of this correlation function shows 
that the parameter c^A corresponds to v^/uqtI where Vh, and r c are, respectively, a sound 
velosity, hard mode frequency at q = and the correlation length of the polarisation in the 
pure lattice. These physical quantities define temperature and concentrational dependence 
of obtained values of a^A. The value of parameter R c /r c (R c is correlation radius of the 
lattice with impurities) defines the ratio of inhomogeneous and homogeneous contributions 
and at T — > T c (R c — > oo ) the line becomes completely inhomogeneous |3|] . In our approach 
the parameter of nonlinearity strongly increases with T — > T c and becomes much greater 
than homogeneous contribution which tends to decrease the magnitude of line maximum. 
Therefore qualititavely our results are in agreement with those obtained in || . 

However comparing the frequency dependence of J(oS) at large uj one can see that it was 
described as (u — f2o)~ 3 / 2 || whereas in our theory as (1 + 4a 2 co , )~ 1 / 2 . More accurate mea- 
surements of line intensity decay could be desirable for both theory comparison. Moreover, 
our theory gives the line shape with two maxima for intermidiate values of a 2 A, which looks 
like that observed in KTN (see fig. 5). Thus the observed line shape with two maxima natu- 
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rally appeares in our theory whereas the second maximum was not obtained in the previous 
theoretical description 0. Its origin was supposed to be some forbidden transition related 
to the mixture of acoustic and optic modes. 

Therefore peculiarities of Raman spectra shape were explained by proposed theory In 
particular it was shown that nonlinear effects leads to strong assymetry of the lines and to 
appearing of the low frequency sharp maximum, its left hand side is defined by dynamic prop- 
erties of the system. The division between dynamic and static characteristics contribution 
into Raman spectra shape makes possible to investigate separately the both aforementioned 
characteristics by Raman spectroscopy method. 
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Figure captions 



Figure 1. Line shape calculated on the base of Eq.(9) for a 2 A =0.22 and (rA) 1= 0.3 
(curve 1), 0.01 (curve 2); doted line is Gaussian form. 

Figure 2. Line shape calculated on the base of Eq.(9) for a 2 A =0.4 and (rA) -1 = 0.1 
(curve 1), 0.05 (curve 2), 0.01 (curve 3); doted line is Gaussian form. 

Figure 3. FOR scattering line shape, solid line - theory at ct 2 A =0.38 and (tA)~ 1 = 0.06, 
crosses - experimental data for KLT with 1% of Li at T = 10 K [3]. Intensity is represented 
in arbitrary units. 

Figure 4. FOR scattering line shape, solid line - theory at ct 2 A =0.55 and (rA) _1 = 0.15, 
crosses - experimental data for KLT with 4% of Li at T = 55 K [3]. Intensity is represented 
in arbitrary units. 

Figure 5. FOR scattering line shape of KTN with 15.7% of Nb at T = 150 K [2] (b); 
calculated line shape for a 2 A =0.21, (tA) _1 = 0.004 (a). 

Figure 6. FOR scattering line shape, solid line - theory at a 2 A =0.25 and (tA) _1 = 0.0036, 
doted line - for KTN with 15.7% of Nb at T = 142 K [2]. 
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